
SUBSTITUTE SPECIFICATION 

Description 

Method and electronic device used to svnthesise the sound of church organ 
flue pipes, by taking advantage f the physical modeling technique of acoustic 
instruments. 

The present patent application refers to a method ad electronic device used 
to synthesise the sound of church organ flue pipes, by taking advantage of 
the physical mod elli ng modeling technique of acoustic instruments. 
Numerous numerical algorithms of physical-mathematical models have been 
developed based o the examination of the physical behavior of organ flue 
pipes and the sound they produce, in order to synthesise the sound emission 
of aerophone instruments in real time. Some of these models are based on 
the mutual symbiotic interaction between a non-linear active section, 
generally defined as "excitation", and a linear passive section, generally 
defined as "resonator". An e xamp le can b e found w i th i n th e m e thod 
d e scr i b e d — in — US — pat e nt — 5,521,328. The relative numerical algorithm 
extemporarily produces a sequence that represents the sound of the 
instrument analysed and translated into a physical model. The sound is 
characterized by an initial time interval, defined as "attack transient", during 
which intensity increases up to a certain value. The intensity value is 
indefinitely maintained over time during the second phase, defined as 
"sustain phase", during which the waveform is approximately periodic. The 
analytical characteristics of this waveform, of which the most important is 
fundamental frequency, depend on each of the parameters that regulate the 
operation of the numerical simulation. Being the simulation performed in the 
time domain instead of the frequency domain because of the presence of 
numerous non-linear functional blocks, the relation between the set of 
parameters and each spectral characteristic of the generated sequence is 
extremely difficult to establish a pirori. 

The characteristics can be altered by charging the set of parameters, 



often empirically, and then evaluating the effect of such a change a 
posteriori. In particular, the fundamental frequency also depends on the 
quantitative characteristics of excitation, and not only on the frequency 
response of the resonator; being the evolution of the sequence extremely 
chaotic during the attack transient phase, the phase of the fundamental 
frequency cannot be pre-determined once the sustain phase has been 
reached. These two peculiarities are unacceptable in high-polyphony 
electronic musical instruments, such as church organs. 

• Other phys i cal - mathemat i ca l mode l 3, 03 the invention dcocribcd in UG 
p at e nt 5,587,548, aro basod on tho conjunct uoago of PCM aud i o 

• synthes i s and physical — ■ math e matical simu l at i on of parts of th o.. 

■ instrument to be reproduced. By ana l ytical l y decomposing th e sound 
camploc of th e instrumontc to bo imitated (or of parts of thorn, ao tho only 
resonant body), and d i viding what can bo eas i ly and cheaply oimu l ated 

-f rom what i s moro conven i ont to storo as part of a wavotablo, a good 

■ compromi eo — b e tw oo n — m e mory — u s ag e — and — computational — powor- 

■ necessary to implement such method can be obta i ned . The excitation 

■ se quence, wh i ch is proproc o ssed by tho algor i thm which simulates port - 
of the acoust i cal bchavoiur of tho instrument (previously analyzed and 

■ math e matica ll y i nterpr e t e d), i s usua ll y storod as a wavotab l o. Howovor, 
■ said method, though requiring computationa l powor for tho physica l 

mathemat i ca l s i mulat i on, implies to s amp l e, ana l yze, and pro - comput e 
th e sound of oach i nstrument to bo reproduced, and said i nstrumen t 's 
-reproduced s ound i s in a ny case bound to said operations, and in 
particular to what is stored i n tho wavotab l oo. 

The present invention consists in an audio-digital synthesis system based 
on digital signal processors, which contains a programme of physical 
simulation of the sound generation of organ flue pipes. The programme is 
divided into three fundamental, conceptually independent sections: the 
first section generates the harmonic part of the sound; the second section 
generates the aleatory part of the sound; the third section processes 
these components by means of a transfer function with two inputs and 



one output, thus obtaining the sequence that represents the sound of the 
organ pipe. Because of the independence of the section that generates 
the harmonic part of the sound, the fundamental frequency and the phase 
of the whole waveform generated by the programme can be determined a 
5 priori. 

The numerical parameters of the simulation programme are partially 
contained in a static memory and partially obtained by processing 
information from an electronic musical keyboard and from a set of user 
controls in real time. They determine the fundamental characteristics of 
10 the generated sound, among which the main characteristics are pitch, 
intensity, time envelope, harmonic composition and aleatory component. 
■ Bo i ng not any information dorivod from rool mus i cal inotrumonto' ooundo 
• and stored as wavotob l oo, memory usage is quito rootrainod. 
For major clarity the description of the method and device according to 
15 the present invention continues with reference to the enclosed drawings, 
which are intended for purposes of illustration only and not in a limiting 
sense, whereby: 

- Figure 1 shows a realisation of a digital electronic musical instrument 
used to synthesise sounds of musical instruments by taking advantage of 

20 the physical modol l ing modeling technique of the invention. 

- Figure 2 shows the three fundamental functional blocks and relative 
interconnections of an audio digital synthesis programme of the sounds 
of church organ flue pipes according to the invention. 

- Figure 3 shows a flow chart that explains one of the three blocks of Fig. 
25 2, according to which a sequence that represents the harmonic part of 

the sounds of church organ flue pipes according to the invention is 
generated. 

- Figure 4 shows a stable realisation of a digital harmonic oscillator with 
two status variables according to the invention. 

30 - Figure 5 shows a procedure used to generate the time variation of the 
operational frequency of the harmonic oscillator shown in Fig. 4 
according to the invention. 



- Figure 6 shows a flow chart used to generate the aleatory component of 
the time progression of the operational frequency of the harmonic 
oscillator shown in Fig. 4 according to the invention. 

- Figure 7 shows an example of time envelope used in the generation of 
5 the sequence that represents the harmonic part of the sounds of flue 

pipes according to the invention. 

- Figure 8 shows a flow chart of a low frequency oscillator used in the 
generation of the sequence that represents the harmonic part of the 
sounds of flue pipes according to the invention. 

10 - Figure 9 shows a time progression composed of non-rectilinear 
sections, according to which the frequency of an oscillator can be 
changed without perceiving an alteration of timbre pitch according to the 
invention. 

- Figure 10 shows an algorithm for the generation of a pseudoimpulsive 
15 periodic sequence according to the invention. 

- Figure 1 1 shows a set of interconnected functional blocks that explains 
one of the three blocks of Fig. 2, according to which a sequence that 
represents the aleatory part of the sounds of church organ flue pipes 
according to the invention is generated. 

20 - Figure 12 shows a status device used to limit the difference between 
two consecutive samples of a sequence according to the invention. 

- Figures 13 and 16 show an example of wave envelope used during the 
attack transient phase of the generation of sounds of flue pipes according 
to the invention. 

25 - Figure 14 shows a wave envelope used to generate the aleatory 
component of the sounds of flue pipes according to the invention. 

- Figure 1 5 shows an architecture that explains one of the three blocks of 
Fig. 2, representing a mathematical model of the resonator of the church 
organ flue pipes according to the invention. 

30 - Figure 17 shows the mutual interaction between two functional blocks 
necessary for the realisation of a' generic harmonic oscillator according to 
the invention. 



- Figure 18 shows an example of a pseudoimpulsive periodic waveform 
generated by the algorithm of Fig. 10, used to generate the aleatory 
component of the sounds of flue pipes according to the invention. 

- Figure 19 explains he operation of the status machine of Fig. 12 according 
to the invention. 

With reference to the aforementioned figures, the electronic musical 
instrument of the invention is physically composed of a set of components, 
whose type, arrangement and interconnection are shown in Fig. 1. 
The embodiment is shown for mere illustrative purposes, since it neither 
represents the central innovation element of the present patent nor the only 
and necessary realization of an electronic musical instrument used to 
synthesise the sound the organ pipes by means of algorithms of physical- 
mathematical simulation. With reference to Fig. 1, the information from a 
musical keyboard (1) and a set of user controls (2) is processed by a control 
unit (3), which regulates the operation of a DSP (6) by means of a plurality of 
numerical parameters contained in a ROM (4). The DSP (6) executes the 
synthesis programme of the sound of the organ pipe in real time, upon 
management from the control unit (3), using a RAM (5) to write and read 
temporary data. The product of the synthesis programme is a numerical 
sequence that is suitably converted by a DAC (7) into the analogue signal 
representing the sound of the organ pipe, which can e reproduced with an 
amplification system and a loudspeaker (8). The synthesis programme, 
which is the central innovation element of the present patent, includes three 
sections. Each section has a fundamental function in the numerical 
simulation of the sound emission of the organ pipe, as shown in Fig. 2. 
The block (9) generates a main harmonic sequence (10) composed of a 
6e ri e s set of harmonic lines, whose amplitude and frequency conveniently 
change over time. By using this sequence and taking advantage of part of 
the rts composition, the block (11) generates a pseudoaleatory signal that 
represents the chaotic component of the sound. The aforementioned 



sequences are the two input signals of the linear resonator (12) that 
models the frequency response of the resonant part of the multiple 
qualities of organ flue pipes, and whose output (13) is the sequence that 
represents the sound of the organ pipes. 
5 The block diagram of Fig. 3 is a detailed view of the functional blocks of 
the harmonic component generator (9). The oscillator (14) generates an 
approximately sinusoidal waveform (16). The fundamental frequency of 
the waveform changes over time within a range of values comprising the 
fundamental frequency of the generated musical note. The details of the 
10 embodiment of the oscillator and the criterion used to change frequency 
over time are illustrated below. 

The waveform (17) is obtained from the sequence (16) through the non- 
linear block (15): if the sequence (16) were exactly a sinusoidal sequence 
x[n] = sin[co 0 n], 

15 the sequence (17) would be 

y[n] = 2 (sin[<B 0 n]) 2 - 1 = -cos[2co 0 n] = sin[2co 0 n - rc/2], 
that is to say a sinusoid with double frequency than the sequence (16). 
Each of the two sequences (16) and (17) is amplified by the relevant 
multipliers (18a) and (18b), and limited by the functional blocks (19a) and 

20 (19b) to values within the +CLIP1 and +CLIP2 intervals. The outputs of 
the blocks (19a) and (19b) are multiplied by two sequences produced by 
the envelope generators (20a) and (20b), respectively, as illustrated 
below, and the resulting products are summed to the node (21). The sum 
is a sequence produced by a series of linear and non-linear 

25 instantaneous operations performed on the waveform (16). If the 
waveform were exactly the sequence x[n], a sequence would be obtained 
in the node (21) whose spectrum would be formed • by harmonic 
components multiple of co 0 (including a> 0 )- 

As illustrated below, the sequence (23) is a low frequency waveform, 
30 whose purpose is the amplitude modulation of the harmonic sequence 
through the product (22). 

The element (24) is a delay line whose impulse response is the sequence 



5.1 [n - N]. Together with the products and the sum of the block (25), this 
element forms a linear filter whose impulse response is 
CBYP + CDELz N . 

The block (26) is a non-linear instantaneous function described by the 
5 following formula: 

f{x) = (x + x 0 ) - (x + x 0 ) 3 + y 0 
where x 0 and y 0 are independent parameters. The purpose of the block is 
to modify the mutual proportion between the amplitudes of the harmonic 
components of the sequence processed by the block. 

10 The block (27) is a band-pass filter, whose peak frequency corresponds 
to the fundamental frequency of the input sequence. The parameter Q of 
the filter is tuned up to obtain the only fundamental frequency of the input 
harmonic sequence with excellent approximation. Moreover, being the 
phase response of the filter null in correspondence of the peak frequency, 

15 the phases of the fundamental frequency of the input and output signals 
of the filter are equal. This characteristic enables to sum the input and 
output sequences of the filter, with no elision effect in the fundamental 
frequency: the block (28) sums the sequences (weighing them with the 
parameters GAIND and GAINF), in order to alter the proportion in 

20 amplitude between the fundamental harmonic component and the group 
of all other harmonic components. The output of the block (28) is the 
main harmonic sequence (10). 

The sinusoidal oscillator (14) consists in a special embodiment of the 
ordinary harmonic oscillator with two status variables, with necessary 
25 measures to improve the robustness to the variation of the operational 
frequency in real time. 

Fig. 4 shows the cycle of operations performed at each sampling interval 
on the two conveniently initialised status variables VAR1 and VAR2. The 
parameter F determines the frequency of the sinusoid produced by the 
30 status variables oscillator that is composed of the steps (29) and (31) in 
the ordinary configuration/The disadvantage of the ordinary configuration 
is that it cannot suffer variations of the parameter F in real time without 



altering the amplitude of the sinusoids described by the same variables, 
in function of the current value of the status variables. Moreover, 
depending on the numerical precision of the oscillator's status variables, 
reductions of the oscillation amplitude can occur even in stationary 
conditions. It is sufficient to amplify the variable VAR2 by a factor 1+e 
(with e positive, but close to zero) by means of the step (30) and limit the 
width of the variable VAR1 by means of the step (32) to values within the 
interval ±1. Using these measures, the variable VAR1 describes a unitary 
amplitude sinusoid with excellent approximation. This variable is the 
output (16) of the block (14) of Fig. 3. The parameter F depends on the 
frequency f according to the relation 

F(f) = 2sin(7if/f sr ) 

where f sr is the sampling frequency. The frequency f can vary in real time 
within an interval [f 0 - Af, f 0 + Af] sufficient to have the frequency changes 
perceived, without a collateral amplitude alteration. 

Having defined the deviation from the central frequency f 0 as 8f, this 
parameter changes in real time according to the scheme of Fig. 5. 
Likewise the signal (23), the signal (33) is a low frequency waveform 
whose purpose is the frequency modulation of the generated sinusoid; 
with the support of the variable VAR1, the block (34) generates an 
aleatory waveform of "sample and hold" type, according to the scheme of 
Fig. 6. Ultimately, 5f varies according to a constant PITCH parameter 
(which, assuming a value in an arbitrary interval [1-5, 1+5], determines 
the fine tuning of the sinusoid) of an oscillating sequence (33) and the 
aleatory sequence (34). The block (34) is described in Fig. 6: every time 
the variable VAR1 passes from a negative value to a positive value, the 
variable RNDPTCH is updated to a new value NEWRND, which is an 
aleatory variable with a probability density function uniformly distributed in 
the interval [1-5RNDP, 1+5RNDP], being 5RNDP an independent 
parameter. 

The two generators (20A), (20B) produce two 5-segment envelope 



signals, whose progression is generically illustrated in Fig. 7. T1...T4 are 
the time intervals in which the signal passes from level LO to L1, from L1 
to L2, from L2 to L3 and from L3 to zero, respectively. The generators 
start producing the respective envelope signals upon a "note on" event. 
5 Level L2 is maintained over time for an indefinite interval SUSTAIN, 
whose end coincides with the corresponding "note off' event. Each of the 
two generators uses its own set of these 8 parameters. 
The signals (23) and (33) are produced by a "Low Frequency Oscillator" 
shown in Fig. 8. The generation method of the triangular waveform with 

10 unitary amplitude and frequency TRFREQ illustrated in the block (35) is 
implicit. The parameters TRFREQ, TRAMPL, TROFFSET, TRCOEFF1 
and TRCOEFF2 determine the conformation of the two signals (23) and 
(33), whose common fundamental frequency is TRFREQ. In particular, 
the signal (32) is a triangular wave of average value TROFFSET and 

15 semi-amplitude TRAMPL, while the signal (33) is formed by sections of 
parabolas, as shown in Fig. 9. The relation between the values 
TRCOEFF1, TRCOEFF2 and the independent parameter K is biunique. 
The special progression of the signal (33) is necessary to obtain a 
triangular frequency modulation as exactly as possible (ref. Fig. 5) around 

20 the nominal frequency f 0f with equal progressions of the positive and 
negative semi-periods, if they are expressed in semitone cents. 
The architecture of the generator (1 1 ) of Fig. 2 is illustrated with details in 
Figures 10, 11 and 12. With reference to Figs. 3 and 10, the signal (16) 
produced by the sinusoidal oscillator (14) is amplified by a factor 

25 RTINGAIN, limited in amplitude by the block (36) to values within the 
interval ±1 , and then processed by the high-pass filter (37). Finally the 
non-linear block (38) cuts the signal's negative values. At the output of 
the block (38) the signal (illustrated in Fig. 13) produced by the envelope 
generator (39) is summed and the result is multiplied by the parameter 

30 RTGAIN. The result RATE is a sequence of values used in the non-linear 
block (42) defined as "RATE LIMITER", which is part of the structure 
described in Fig. 11. With reference to Fig. 11, the functional block (40) 



generates a white aleatory sequence, with a uniformly distributed 
probability density function processed by the low-pass filter (41). The 
obtained sequence is the input signal of the structure formed by the delay 
lines NBDL1, NBDL2, NBDL3, NBDL4, the sums NBS1, NBS2, NBS3, 
5 the multipliers NCGAIN, NBFBK and the non-linear block (42). The set 
formed by these elements, including the topology of interconnections, is 
defined as "NOISE BOX". The signal generated by the block (42), which 
is the output of the aforementioned set, is amplified by a factor NGAIN 
and multiplied by the signal produced by the envelope generator (43), 
10 whose time progression is illustrated in Fig. 14. The signal NOISE is the 
output of the generator (1 1 ) of Fig. 2. 

Fig. 12 describes the non-linear block "RATE LI MITER" (42) formed by 
the sums RLS1 , RLS2, the limiter (44) and the unit delay element (45). 
The value memorised in the delay (45) is subtracted by means of the 

15 adder RLS1 from the input signal "IN"; the result is then limited to values 
within the interval ±RATE (being RATE the sequence generated by the 
network illustrated in Fig. 10), and finally summed again to the current 
delay value (45) at the node RLS2. The result "OUT" is memorised in the 
delay element (45) for a successive cycle. Fig. 13 shows the time 

20 progression of the envelope generated by the block (39): upon a "note 
on" event, starting from the level NBLO, the level NBL1 is reached in a 
time NBT, indefinitely sustained over time, also after the corresponding 
"note off" event. Fig. 14 shows the time progression of the sequence 
generated by the block (43): upon a "note on" event, the signal starts 

25 from the value NLO, reaches the value NL1 over a time NT1 and the level 
NL2 over a time NT2 sustained until the successive "note off" event. 
Upon this event the signal reaches the value zero in a time NT3. 
With reference to Fig. 11, the non-linear block "RATE LIMITER" (42) can 
be replaced with a linear filter, whose gain has a progression described 

30 by the same sequence RATE generated by the architecture of Fig. 10, so 
that the structure "NOISE BOX" of Fig. 1 1 is a linear time-variant filter. 
With reference to Fig. 2, the outputs of the generators (9) and (11) are 



the inputs of the resonator (12) illustrated with details in Fig. 15. The 
functional blocks of the network (12) form a cycle of operations, along which a 
sequence of samples propagates for a potentially infinite time. The two 
contributions of the two generators (9) and (11) are added to this sequence, 
5 instant by instant in the sum nodes (46) and (48) nodes, respectively, to 
sustain the energy of the computed sequence. The structure of Figure 15 is 
the translation into a mathematical model of the resonant part of the organ 
flue pipe, defined as "pipework". In particular, the low-pass filter (47) 
emulates the dissipation of acoustic energy, with variable intensities in 

1 0 function of the frequency; the high-pass filter (49) attenuates all the frequency 
components lower than the fundamental frequency; by means of the product 
(51), the envelope generator (50) produces a signal that represents the time 
progression of the loop gain of the resonant system; the filter (52) alters the 
sequence phase, leaving its module unchanged; the factor TFBK (53) 

15 depends on the type of acoustic termination at the top of pipework; finally, the 
delay line BDELAY (54) considers the time needed by an acoustic pressure 
wave to cover the pipework from the base to the top and vice versa. The time 
progression of the signal produced by the envelope generator (50) is traced in 
Fig. 16; likewise the envelope of Fig. 13, upon a "note on" event, the signal 

2 0 passes from a value FBLO to a value FBL1 in a time FBT, and then remains 

constant. The output sequence (13) is the signal emitted by the mathematical 
model of Fig. 2 as a whole, that is to say the time representation of the sound 
emission of the organ flue pipes. 

The description continues with the original innovative characteristics of the 
25 audio digital synthesis technique of the sound of flue pipes. 

The literature on the generation of sounds of instruments with continuous 
sound emission, among which aerophone instruments, by means of the 
physical mod elli ng modeling technique, proposes solutions based on a 
mutual interaction between a non-linear active part, normally defined as 

3 0 excitation (55), and a linear passive part, defined as resonator (56), according 

to the scheme of Fig. 17. Th e m e thod e xpos e d in US pat e nt 



■ 5,521,328 can bo conoidorod as an oxamplo of th i s technique . In the 
case of aerophone instruments, the energy contributed to the system is in 
the form or sound pressure and the signal produced is the progression of 
the sound pressure wave irradiated by one or more suitable points of the 
5 resonator. The waveform p(t) is the progression of the air pressure that 
the performer (or the bellows, in the case of a church organ) exercises on 
the instrument mouthpiece. - According to this progression and to the 
progression of the pressure w(t) in a suitable point inside the resonator, 
an oscillating acoustic pressure e(t) injected in the resonator is 

10 generated. Once the sustain phase has been reached, the pressure e(t) 
has the same fundamental frequency as the pressure w(t). Being linear 
(except for very special operation modes), the resonator can be 
described with an impulse response r(t), which generates the return 
signal w(t) and an impulse response h(t), which generates the output 

15 signal y(t). The latter is the time progression of the sound emission of the 
instrument. Being it a numerical simulation performed in the time domain 
instead of the frequency domain, the fundamental frequency of the 
oscillation on which the system stabilises, once the sustain phase has 
been reached, is extremely difficult to predict mathematically. This 

20 depends on the fact that the frequency depends on the time progression 
of the forcing signal e(t), and not only on the frequency values in which 
the amplitude spectrum of the impulse response of the resonator has the 
relative maximum values. In fact, any type of harmonic oscillator 
(electronic, mechanical, etc.) has this characteristic. With regard to wind 

25 instruments (including organ pipes), it is sufficient, for example, to 
increase the sound pressure to obtain an increase of the fundamental 
frequency of the acoustic wave, in addition to an intensity increase, 
although the characteristics of the resonant part remain unchanged. 
Another inevitable characteristic of the oscillating systems illustrated in 

30 Fig. 17 is the unpredictability of the phase of the generated signal, once 
the sustain phase has been reached. Since the waveform p(t) used to 
stimulate the system is partially chaotic, and in any case it does not 



contain any information about the phase of the stationary wave sustained 
by the resonator, the attack transient of the signal y(t) is always and 
unpredictably different. Therefore, although the waveform has always the 
same periodic time progression in sustain conditions, it is impossible to 
5 determine the evolutions that bring the system towards this progression. 
In quantitative terms, it is impossible to determine the phase of the 
fundamental frequency of any signal processed inside the stable 
oscillating system of Fig. 17, taking the instant when the stimulus p(t) 
starts as time origin. Together with the difficulty encountered in 
10 determining the fundamental frequency a priori, this is unacceptable in 
the field of high-polyphony electronic musical instruments, such as 
church organs. 

The synthesis system used in the present invention derives directly from 
the synthesis in the time domain described in general and is 

15 characterised by the total autonomy of the excitation signal from the 
signal produced by the resonator. In fact, the main harmonic sequence 
(10) extemporarily generated by the block (9) of Fig. 2 is the imitation, as 
faithful as possible, of the signal e(t) of the system of Fig. 17 (assuming 
that the latter is a good mathematical model of the flue pipe of a church 

20 organ), with the substantial difference that the fundamental frequency 
and the phase of this sequence, and consequently of the sequence 
produced by the system (13) as a whole, are perfectly determined a 
priori. 

The preparation of the numerical parameters of any oscillating system, as 
25 generically illustrated in Fig. 17, requires special sensitivity and skill, 
apart from the perfect knowledge of its mathematical model. This means 
that the good operation of the system may be impaired, and the system 
may become unstable or even inharmonious, if only one of the 
parameters has a value not included in a proper range. Moreover, 
30 different operational modes of the oscillator can be obtained only by 
acting simultaneously and with special attention on a plurality of 
parameters, with the risk of making the time evolution of one or more 



signals in transit along the functional blocks of the system uncontrollable. 
This makes the search for multiple sounds produced by this type of synthesis 
slow and difficult. On the contrary, a system with no feedback between 
resonator and excitation, such as the system shown in Fig. 2, enables to 
modify the numerical parameters of the three functional blocks (9), (11), (12) 
in a completely independent way, without impairing the good operation of the 
system as a whole. This allows obtaining a larger variety of sounds than he 
one obtained by means of a feedback loop system with equal complexity. 
Tho current litoraturo proposes, as e.g. in US patent 5,587,5 - 18, an alt e rnat i v e 
t e chniqu e wh i ch i s known as commut e d synth e sis, ba se d on e xc i tat i on 
wav e tab le s and — r e sonant f il t e r s , — b ei ng th e latt e r us e d to s i mulat e th e 
acoust i ca l b e hav i or of an instrum e nt's li n e ar and passiv e part. — I n such a 
cas e , w i th an ad e quat e samp le d sound' s ana l ys i s and optim i zat i on, a - good 
comprom i s e is found b e twe e n th e n e c es sary amount of m e mory for th e 
wav e tab le s — and — the — n e c e ssary — computat i ona l — pow e r — te — impl e m e nt th e 
phy si ca l math e mat i ca l a l gor i thm which corr e sponds to th e instrum e nt's 
r e sonant parts. 

The system of Fig. 3 shows a sequence of operations performed on the signal 
produced by the sinusoidal oscillator (14). The type and order of the 
operations are only one of the possible realizations used to generate a 
waveform sufficiently rich in harmonic components and provided with a 
suitable time evolution. In any case, some of the functional blocks of the 
system, such as the delay (24) and the non-linear function (26), derive from 
mathematical models of wind instruments known in the literature, without the 
need of using them. 

The originality of the system mainly consists in the adaption of an ordinary 
oscillator with status variable to non-stationary operational conditions, by 
developing the functional blocks (30) and (32) of Fig. 4, in order to make the 
oscillator robust to the variations of the parameter F 2 of the block (29). 
With reference to Figures 5 and 6, the originality derives from the 
development of the generator (34) to obtain pleasing random frequency 
variations in real time. Assuming the factor (33) as constant, that is to say 
assuming the absence of the low frequency oscillation of the sequence of, the 
latter assumes a new random value at every period of the sinusoidal 
sequence VAR1 . The result is a statistic uniformly distributed 



variation of the wave period, in terms of probability density function. The 
variation is perceived as a pleasant irregularity in the sound emission. 
Otherwise, if 5f assumes a new random value at every sampling instant, 
the length of every wave period will be described by a variable formed by 
5 the sum of N aleatory contributions, each of them provided with uniformly 
distributed probability density (N is the number of samples per period). In 
view of the Central Limit Theorem, the higher is N, the more the 
probability density function of this variable approaches a Gaussian 
function. The frequency variation would be very irregular, since high 
10 frequency deviations would be obtained much more rarely than small 
deviations from the nominal frequency. This would be very unpleasant, 
since wave periods with very different length from the nominal length 
could be generated and perceived as sudden malfunctions of the 
generation model. 

15 The generator of the aleatory component (11) of Fig. 2 is completely 
original, and the embodiment of Fig. 10, 11 and 12 derives from the 
analysis of samples of sounds emitted by a large variety of organ flue 
pipes, and from some hypotheses on their operation physics. In 
particular, by analysing the spectrogram of the individual wave periods of 

20 a sample and using a much finer time resolution than a wave period, it 
can be noted that a large percentage of sound energy concentrates in a 
time interval much shorter than the period, always situated in the same 
position along the wave period. Such sound energy covers a frequency 
interval considerably higher than the interval covered as an average by a 

25 plurality of periods. Therefore, the characteristics of the spectrogram of 
the stationary part of the sound of flue pipes are similar to the 
spectrogram of a train of equidistant impulses, with the energy of the 
individual period concentrated in each impulse. These considerations 
justify the architecture illustrated in Fig. 10: the sequence RATE is 

30 obtained through a series of elementary deterministic operations 
performed on the sinusoid (16). Once the sustain phase has been 
reached, the sequence RATE, assumes a qualitatively impulsive 



progression, of which Fig. 18 shows one example, where T 0 is the period 
of the sinusoid (16). Regardless of the method used to obtain a 
pseudoimpulsive sequence, the sequence is conceptually one of the 
inventive foundations of the generator (1 1). 

The structure illustrated in Fig. 11 is formed by the four delay lines 
NBDL1, NBDL2, NBDL3, and NBDL4. Together with the sums NBS1 , 
NBS2, and the product NCGAIN, the first three delay lines form a FIR 
filter. The output of this filter (that is to say the sum NBS2) is processed 
by the non-linear element (42) and then, after being multiplied by NBFBK 
and after passing through the fourth delay line and the sum NBS3, 
injected again in the aforementioned filter. If it weren't for the element 
(42), the structure "NOISE BOX" would be a linear filter, whose spectrum 
would have a voluntarily inharmonious progression, with a plurality of 
resonance peaks distributed in a non-deterministic way, depending on 
the length of the delay lines and the two independent parameters 
NCGAIN and NBFBK. These four quantities are dimensioned in order to 
imitate the frequency response of a resonator with irregular geometry, 
such as the portion of space of the organ pipe immediately inside the 
mouth. Because of the periodic oscillation of the sequence RATE, the 
element (42) causes a continuous periodic variation over time of the 
"gain" (not strictly definable as such, since the "RATE LIMITER" is a non- 
linear block) of the entire "NOISE BOX". In particular, with reference to 
Fig. 18, when the sequence RATE assumes the minimum value, the non- 
linear distortions caused by the block (42) imply energy losses that 
heavily reduce the resonance effects of the "NOISE BOX". Vice versa, 
during the (much shorter) instants in which the sequence RATE assumes 
relatively high values, the resonant effect of the "NOISE BOX" emerges 
and the intensity of the aleatory component increases. It can be noted 
that during the attack transient, because of the envelope generator (39), 
whose progression is shown in Fig. 13, the sequence RATE assumes 
higher values than during the sustain phase; this increases the 
resonance of the "NOISE BOX" during the first instants of synthesis, in 



order to simulate the acoustic phenomena defined as chiff, cough, etc. 
produced by the flue pipes if the valve that regulates the passage of air 
from the bellows to the foot is opened rapidly. 

The non linear block (42) is formed by the two adders RLS1, RLS2, the 
limiter (44) and the unit delay element (45). At every sampling instant, the 
difference between the previous output value and the current input value 
is first limited in width to values within the interval +RATE and summed 
again to the previous output value, thus obtainwig the current output 
value. The output sequence "follows" the input sequence, maintaining an 
inclination limited according to the value RATE. For mere illustrative 
purposes, Fig. 19 shows a chart of an input sequence (continuous line) 
and an output sequence (dotted line). In the instant t 0 the inclination of 
the sequence IN exceeds the value RATE/sample, therefore the 
sequence OUT separates until it re-joins at point tj, after which the 
sequence IN remains constant. In the instant t 2 the excessive inclination 
of the sequence IN causes the immediate separation of the sequence 
OUT up to the re-conjunction point U- With respect to a linear filter, the 
advantage of the "RATE LIMITER" is the elimination of possible 
discontinuities of the aleatory sequence, while still maintaining a sufficient 
bandwidth, being such discontinuities extremely unpleasant for the 
human hearing. This aspect represents the originality of the "RATE 
LIMITER". 

The non linear block (42) can be replaced with any functional block 
whose effect on the structure "NOISE BOX" of Fig. 11 is the quantitative 
resonance variation generated by the structure, according to a periodic 
progression. 

As regards the linear resonator (12), the physical considerations that 
involve the choice of the functional blocks of Fig. 15 are described herein. 
The resonant part of an organ pipe, defined as pipework, can be 
mathematically described, in the most elementary way, with a "comb" 
filter 1/(1 - FBK-z" N ), in which the feedback coefficient FBK is related to 
the loop gain of the filter and the parameter N is inversely proportional to 



the first resonance frequency of the same. The more complex resonator 
of Fig. 15 derives from this base, which is very used in the field of audio 
digital processing. Among the elements of the resonator, the function of 
the delay line (54) appears evident. The response in module of the low- 
5 pass filter (47) is designed so as to consider the different energy losses 
suffered by the various harmonic components during their transit along 
the pipework, while the high-pass filter (49), whose cut-off frequency is 
lower than the fundamental frequency of the resonator, completely 
eliminates the continuous component of the stationary wave, to take into 

10 account the fact that the average acoustic pressure inside a pipework is 
approximately equal to the external pressure. Because of the envelope 
generator (50), during the first operation phase of the resonator, the loop 
gain of the system is moderately overabundant with respect to the value 
once the sustain phase has been reached, in order to obtain a faster 

15 initial energy accumulation in the resonator, that is to say a faster attack 
transient of the generated sound. The sign of the factor TFBK (53) is 
especially important: a positive sign for a pipework open at the mouth and 
on top, and a negative sign for a pipework open at the mouth and closed 
on top. This derives from the physics of the reflection of an acoustic 

20 pressure wave in correspondence of the pipework terminations. This 
physical law also justifies the use of the all-pass filter (52), the most 
important element of the resonator from the conceptual point of view. If, 
on one side, the mono-dimensional model of the pipework is sufficiently 
accurate to justify the use of an individual delay line to simulate the 

25 longitudinal propagation of an acoustic wave in the pipework, the 
approximation becomes unacceptable in the wave reflection in 
correspondence of a structural discontinuity characterized by non- 
negligible transversal dimensions, such as the top of the pipework. The 
all-pass filter (52) modifies the total phase delay of the closed cycle 

30 formed by the elements (46) ... (54) in a selective way with respect to the 
frequency, in order to make the' resonance of the linear resonator (12) 
realistically inharmonious. The same filter is optionally used to modify the 



value of the first resonance frequency of the pipework in real time through 
controlled variations of its coefficients upon a "note off' event, to simulate 
the phenomenon of the moderate reduction of the fundamental frequency 
of the sound of small flue pipes when the air inlet valve closes. 



